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Abstract 

We calculate the Emptiness Formation Probability (EFP) in the spin-Calogero 
Model (sCM) and Haldane-Shastry Model (HSM) using their hydrodynamic de- 
scription. The EFP is the probability that a region of space is completely void 
of particles in the ground state of a quantum many body system. We calcu- 
late this probability in an instanton approach, by considering the more general 
problem of an arbitrary depletion of particles (DFP). In the limit of large size 
of depletion region the probability is dominated by a classical configuration in 
imaginary time that satisfies a set of boundary conditions and the action calcu- 
lated on such solution gives the EFP/DFP with exponential accuracy. We show 
that the calculation for sCM can be elegantly performed by representing the 
gradientless hydrodynamics of spin particles as a sum of two spin-less Calogero 
collective field theories in auxiliary variables. Interestingly, the result we find for 
the EFP can be casted in a form reminiscing of spin-charge separation, which 
should be violated for a non-Hnear effect such as this. We also highhght the 
connections between sCM, HSM and \ = 2 spin-less Calogero model from a 
EFP/DFP perspective. 
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1. Introduction 

One-dimensional integrable models have an important role in the study of 
strongly correlated systems. When the reduced dimensionality makes interac- 
tion unavoidable, perturbative techniques can quickly loose applicability and 
over the years more sophisticated tools have been developed to tackle these 
problems. These tools clearly involve certain approximations and the existence 
of an exact solution for some models can allow to check their validity. 

The conventional approach in solving quantum integrable model is known 
as Bethe Ansatz (and its generalization). It is very successful in constructing 
the thermodynamics of a system, but not very suitable to study its dynamics 
and the correlation functions, due to the increasing complexity of its solutions. 
However, a very elegant formalism was developed using the Quantum Inverse 
Scattering Method (QISM) [H to express correlation functions as determinants 
of certain integral operators (Fredholm determinants). In this formalism, the 
simplest correlation function one can write is known as the Emptiness Formation 
Probability (EFP) and measures the probability P{R) that a region of length 
2R is completely void of particles. For lattice models, one is interested in P{n), 
the probability that n consecutive lattice sites are empty. In spin chain, taking 
advantage of the Jordan-Wigner mapping between particles and spins, the same 
quantity can be thought of as the Probability of Formation of Ferromagnetic 
Strings (PFFS), i.e. the probability that n consecutive spins are aligned in the 
same direction. 
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One should notice that the EFP is an n-point correlator and is, therefore, 
a much more complicated object compared to the usual two-point correlation 
functions one normally studies in condensed matter physics. However, due to 
the strongly interacting nature of the 1-D model, the QISM tells us that it is 
in fact no worse than other correlators between two points a length n apart 
and even somewhat simpler and more natural. Moreover, the EFP is one of 
those extended objects like the Von Neumann Entropy, or the Renyi Entropy, 
that in recent years have attracted a lot of interest because of their ability to 
capture global properties that were not observed before from the study of 2-point 
correlation functions. The latter quantities are of course motivated by studies 
of entanglement and quantum computation, while the EFP arises naturally in 
the contest of integrable theories. 

Despite the claimed simplicity, the calculation of the EFP is by no means an 
easy task. For some models, the specific structure of the solution has allowed to 
find the asymptotic behavior of the EFP as n ^ oo. For instancejthe EFP in the 
whole of the phase-diagram of the XY model was calculated in [3, [3, 3 01 using 
the theory of Toeplitz determinants, while for the critical phase of the XX Z spin 
chain the solution was found in 0,0,131 using a multiple- integral representation. 
The EFP has been considered also for the 6-vertex model Hi, for higher 

spins XXZ and for dimer models [U]. We also remark that high temperature 
expansions of the EFP for Heisenberg chains have been studied in 
A recent review of the EFP can be found in (H| or [l3| . 

Field theory approaches are normally most suited for the calculation of large 
distance asymptotics of correlation functions, but conventional techniques like 
those inspired by the Luttinger Liquid paradigm (i.e. bosonization) are not 
appropriate for extended objects like the EFP and only capture its qualitative 
behavior, while being quantitatively unreliable, as it was showed in The 
reason for this failure is that Luttinger Liquid is applicable only to low-energy 
excitations around the Fermi points, where the linear spectrum approximation 
is valid, while correlators like the EFP involve degrees of freedom very deep in 
the Fermi sea, where the whole spectrum with its curvature is important. 

For this reason, the field theory calculation of the EFP requires a non-linear 
generalization of conventional bosonization, i.e. a true hydrodynamic descrip- 
tion of the system. In [r?| it was shown that, with such a non-linear collective 
description available, the calculation of the EFP is possible by employing, for 
instance, an instanton approach. 

In this paper, we will extend the machinery developed in [l3| and apply 
it to the spin-Calogero Model (sCM), for which a (gradientless) hydrodynamic 
description was recently constructed from its Bethe Ansatz solution [19||. The 
sCM is the spin— 1/2 generalization fi^, [21I, of the well-known Calogero- 
Sutherland model and is defined by the Hamiltonian 



where Pj-; is the operator that exchanges the positions of particles j and I. We 
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chose to analyze this Hamiltonian assuming it acts on fermionic particles, which 
means that the exchange term selects an anti-ferromagnetic ground state (l9| . 
The coupling parameter A is taken to be positive and N is the total number of 
particles. 

In , a collective description of the model was derived using four hydrody- 
namic fields: the density of particle with spin up/down p^^ i and their velocities 
U|,X. The Hamiltonian in terms of these fields is valid only for slowly evolv- 
ing configurations, where terms with derivatives of the density fields can be 
neglected. This description is referred to as a gradientless hydrodynamics. In 
[l9j , this theory was used to show the non-Hnear coupling between the spin and 
charge degrees of freedom beyond the Luttinger Liquid paradigm and it was 
shown that, while a charge excitation can evolve without affecting the spin sec- 
tor (for instance for a spin singlet configuration), a spin excitation carries also 
some charge with it, in a non-trivial way. 

The EFP for the sCM has not been considered in the literature yet. For the 
spin-less case of the Calogero-Sutherland interaction, the asymptotic behavior 
of the EFP was obtained using the form of the ground state wavefunction and 
thermodynamical arguments [i^l (see (1, [l3| for details). It should be noted 
that for certain special values of the coupling parameter A, the spin-less theory 
is tightly Hnked with Random Matrix Theory (RMT) and the EFP is the prob- 
ability of having no energy eigenvalues in a given interval. For these values of 
A the EFP can be calculated with much greater accuracy due to the additional 
structure provided by RMT 

If we write the ground state of the system as '^g{xi,X2, . . . ,xn), the Empti- 
ness Formation Probability is defined as 

PiR) = L N / dx,...dxN\^G{xi,...,XN)\\ (2) 

\^G\^G/ J\xj\>R 

or, following [H 

P(i?)= lim (*G|e""^"«''^^"^'^"|*G) , (3) 

a — >oo 

where pd'x) is the total particle density operator 

N 

Pc{x) = ^S{x-Xj) . (4) 
i=i 

For a model like the sCM, we can also introduce the EFPs for particles with 
spin up or down separately 

P^.i(R) = lim (*G|e~"^-"«''T-i(")^"|*G) , (5) 

which will allow us to discuss the EFP as well as the PFFS. 

The approach we use to calculate the EFPs ^ in this work is similar to what 



was explained in [17(1 . The idea is to consider the system as a quantum fiuid 
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evolving in imaginary time (Euclidean space) . Then the EFP can be considered 
as the probability of a rare fluctuation that will deplete the region —R < x < R 
of particle at a given imaginary time (say t = 0) . With exponential accuracy, 
the leading contribution to this probability comes from the action calculated on 
the saddle point solution (instanton) satisfying the EFP boundary condition. 

In section [2] we will flrst review the results of and transform them into 
an intriguing form where the dynamics can be decoupled into two independent 
fluids of spin-less Calogero-Sutherland particles. This two-fluid description is 
one of the interesting observations of this paper. In section [3] we will explain 
the instanton approach and formulate the problem in this language. In section 
[4] we will concentrate on a generalization of the EFP, the Depletion Formation 
Probability (DFP) which was introduced in This correlator will allow us 
to calculate the different EFPs very efficiently by taking its different limits in 
section [H Most noticeably, we will derive the PFFS for the Haldane-Shastry 
model as the freezing limit of the sCM. In section [6] and [7| we will consider two 
additional DFP problems. Instead of specifying boundary conditions for both 
the spin and charge sectors of the fluid as we did in the previous sections, we will 
now relax these conditions and constrain only one component at a time: this 
analysis suggests that an effective spin-charge separation can be conjectured for 
the EFP/DFP of the sCM. In section [8] we combine all these results and suggest 
a physical interpretation of them. The flnal section contains some concluding 
remarks. To avoid interruptions in the exposition, certain technical formalities 
are moved to the appendices and are organized as follows. In appendix [Al we 
will revise and adapt the calculation of [1/1] to calculate the instanton action 
for our cases. In appendix [B] we will repeat this calculation in the linearized 
hydrodynamics approximation or bosonization, to aid the discussions in section 



2. Two-fluid description 

In jlil the gradientless hydrodynamic description for the sCM ^ was de- 
rived in terms of densities and velocities of spin up and down particles: p-\,i{t, x), 
v^^lit, x). Here, we prefer to use densities and velocities of the majority and mi- 
nority spin: pi,2(i, a:^), ^'i,2(i, a;), i.e. the subscript 1 (2) takes the value t or | 
which ever is most (least) abundant species: 

_ PT +Pi + Ipt -Pil Pc + Ps 
PI = ^ = ' 

_ PT +Pi ~ \P^ - Pi\ Pc - Ps ,„s 

where we introduced the charge and spin density 

Pcit.x) = p^+pi=pi+p2, (8) 
Psi.t,x) = Ip-^ - pi\ = pi - P2 . (9) 
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Please note that whatever species is majority or minority is decided dynamically 
in each point in space and time. 
Under the condition IigII 



\Vi -V2\< TTPs 



the Hamiltonian is 
1 



H 



127r(A + 1) 



-foo 



dx 



"■LI 



2A + 1 



1,3 ^ 



where 



kRi,Li 

kR2,L2 



= Vi±{X + l)7rpi ± \ttp2 , 

= (A + l)v2 - Xvi ± (2A + l)7rp2 



(10) 



(11) 



(12) 



are the four dressed Fermi momenta. 

It turns out that an auxiliary set of hydrodynamic variables decouples the 
Hamiltonian (fTTj) into the sum of two independent spin- less Calogero-Sutherland 
fluids a and b: 







Ax 



where 



Pa 



Pb 



Vb 

Xa 
Xb 



a—a,b '' 



kRl - k Ll 

2TTXa 
fcj?.2 - fc£2 

27rAfc 
kRl + kLl 



1 



2 \ 2 



TT^A 



-pI 



X 



= Pi 



1 



A + 1 

P2 , 



P2 



«J?,2 



«L2 



A+1 
(A + 1) V2 



Xvi 



A + 1 , 

(A + 1)(2A+1) 



(13) 

(14) 
(15) 
(16) 

(17) 

(18) 
(19) 



We remark that both the auxiliary variables and the real variables satisfy 
the canonical commutation relations, i.e. 



[pa{x),vi3{y)] = -ih6a,i3S'{x -y) , a, (3 ^ {1, 2}; {a, b} 



(20) 



The form of the Hamiltonian (flSl) is one of the interesting observation of this 
paper, since it allows us to reduce the spin Calogero-Sutherland model into a 
sum of two spin- less theories. Each of the terms in square brackets in (fTSl is the 
gradientless Hamiltonian of a spin-less CS system with coupHng constants Xa.b 
given by (fTSl fT9|) . In [l3| the gradientless hydrodynamics of spin-less particles, 
like the ones in ifTSj) was used to calculate the EFP from the asymptotics of an 
instanton solution. In the next section we review this approach and we leave 
the mathematical details to appendix [Al 
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3. The instantonic action 



Let us perform a Wick's rotation to work in imaginary time r = it. Note 
that this makes the velocities in (fT2|) imaginary {v iv) and the fc's complex 
numbers. The x — t plane is mapped into the complex plane spanned by z = 

X + IT. 

Following [l3|, we will calculate the EFP as an instanton configuration (i.e. 
a classical solution in EucHdean space) that satisfies the boundary condition 

Pa{T = 0;~R< X < R) = , a=l,2,c, (21) 

in the limit i? ^ oo, i.e. R much bigger than any other length scale in the 
system. This hmit guarantees that the gradient-less hydrodynamics ifTS]) is 
valid in the bulk of the space-time. Once we have the classical solution of the 
equation of motion 0efp that satisfies l(2T|) . a saddle-point calculation gives 
the EFP with exponential accuracy as the action S calculated on this optimal 
configuration \v^ : 

P{R) ~ c-^i'i'^'"'^ . (22) 

Of course, to uniquely specify the problem, the boundary conditions at infinity 
have to be provided as well and we will take them to be those of an equilibrium 
configuration: 



Pi{t,X) poi , Vl[T,X) 0, 

P2[t,x) -> po2 , V2[t,x) 0. (23) 



When poi = po2 we have an asymptotic singlet state (the AFM in zero magnetic 
field) . The condition poi 7^ P02 can be achieved via a constant external magnetic 
field which would result in a finite equilibrium magnetization. It is easy to 
implement these boundary conditions in our two-fluid description using (|14m7p . 

The key point for the calculation is that we can represent the hydrodynamic 
fields in terms of the dressed Fermi momenta km, fcfl2 (which in Euclidean space 
become complex and complex conjugated to fcLi, ^L2 respectively) through (fT2|) : 

km = Xairpa + iVa , kB.2 = AfcTTpf, -|~ \Vb . (24) 

In [l^, it was shown that these fc-fields propagate independently according to 
4 decoupled Riemann-Hopf equations 

drW - iwdxW = , w = kR^L-1,2 ■ (25) 

These equations have the general (imphcit) solution 

w = F{x + iwt) (26) 

where F{z) is an analytic function to be chosen to satisfy the boundary condi- 
tions. 



7 



Guided by [17||, the solution for an EFP problem is 

km = Fa{x + ifc^ir) , kR2 Fb{x + ifcfl2T) , (27) 



with 



Fb{z) = XbTTpm + Afc7r?7b , ^ 1 ' (29) 

V V — it^ / 

which automatically satisfy the conditions at infinity (|23|) : 

Pa(T, a; -> oo) poi + T— — rP02 = POa , 

A + i 

P6(r, a; -> cx)) ^71-^^02 , = Pob , 

V a, b{T ., X Co) 0, (30) 

while 'qa,b are two, possibly complex, constants that allow to satisfy the EFP 
boundary conditions IpTI) . 

In appendix [A] we show that the instanton action can be expressed as a 
contour integral where only the behaviors of the solutions (|27| at infinity and 
close to the depletion region are needed, saving us the complication of solving 
the implicit equations in generality. Using the two-fiuid description, the action 
can be written as the sum of two spin-less Calogero-Sutherland fluids: from l(96l) 
we have ^ 

5efp = 2 '^^^^ X! ^" ^" ■ (^-"^^ 

Oi—a.h 

Before we proceed further, we should mention that the two-fluid descrip- 
tion we employ is valid as long as the inequality (fTO| is satisfled. In fact, the 
solution l(27|) could violate the inequality in a small region around the points 
(r, x) = (0, ±i?). However, close to these points the hydrodynamic description is 
expected to be somewhat pathological, because gradient corrections (which we 
neglect) become important. As it was argued in [3, |T3|, the contributions that 
would come to the EFP from these small regions are subleading and negligible, 
in the asymptotic limit i? ^ 00 we consider. Therefore, we do not need to worry 
about what happens near the points (r, x) = (0,±_R). However, a consequence 
of the "singular" nature of these points is that, in our solution, the species that 
constitutes the majority (minority) spin in the region of depletion —R < x < R 
at r = 0, could switch and become minority (majority) at inflnity. This could 
be important to keep in mind in interpreting our formulae, but our formalism 
already takes that into account naturally. 

4. Depletion Formation Probability 

It is more convenient to consider a generalization of the EFP problem, called 
Depletion Formation Probability (DFP) which was introduced in In hydro- 
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dynamic language the DFP boundary conditions for the majority and minority 
spins are 

Pi{t = 0;-R< X < R) 



P2{t = 0; -R < X < R) 



Pi 

P2 



(32) 



The DFP is a natural generalization of the EPF l(2T|) and it reduces to it for 
Pi, 2 = 0. Of course, there is some ambiguity on the microscopic definition of 



the DFP (see 



One can, for instance, consider it as the macroscopic 



version of the s-EFP introduced in [25] . We will first calculate the DFP as the 
most general case and later take the appropriate interesting limits. 

In terms of the auxiliary fields we introduced in (|14m7p to achieve the two- 
fiuids description l|13p . the DFP boundary conditions are 

Pair = 0, -R < X < R) = pi + -r/°2 = Pa 



Pb{T = 0, -R < X < R) 



1 



A + 1 



A + 1' 



P2^ Pb ■ 



(33) 



We specify the parameters 7/1,2 in l|28|29p by expressing the boundary con- 
ditions l|33p in terms of the dressed momenta using 

A 



Pi 

P2 

Vl 
V2 



1 

7r(A + l) 
1 

7r(2A + 1) 
Im km , 
1 



Re kfii 
Re A:_R2 



2A+ 1 



Re kR2 



A + 1 



Im kji2 + A Im km 



(34) 



This leads to 



Aq77q 
Ab?76 



Aa (POa - Pa) 

(A + l)(poi-pi) + 
Ab {pob - Pb) 

(2A+I)(p02-P2) 



Hpo 



P2) 



(35) 



It is now straightforward to obtain the DFP by substituting l(35|) into ([ST 
After some simple algebra we get 



DFP 



(R) 



exp< 



2 r 



A [poc 



Pcf + (POI 



Plf + iP02 - P2Y 



(A + 2) (Poc 



Pcf + 2 (POs 



~ n2 
Ps) 



(36) 



where we introduced a notation in terms of the charge field (poc = Poi + P02, 
Pc = h + P2) and of the spin field (pos = Poi - P02, Ps = Pi ^ h)- Equation 
l(36|) is the main result of this work. To understand it better, we will consider 
several interesting limits. 
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4.1- Asymptotic singlet state 

If no external magnetic field is applied, the equilibrium configuration of an 
anti-ferromagnetic system like the one we consider is in a singlet state. This 
means that in the boundary conditions at infinity l(23l) we should set poi = P02 = 
po- In this limit l(36|l reduces to 



nSinglet . I ""^ 

Pdfp (R) = "^""P] Y 



A {2po - pcf + (po - Pif + (po - hf 



exp-^ 



TT 



{X+l){2po-pcf + ^pl 



R' 



(37) 



5. Emptiness Formation Probability 



By taking the limit pi^2 = we can use ((36l) to calculate the different EFPs. 
The probability to find the region —R <x<RatT — Q completely empty of 
particles is therefore 



EFP 



[R] 



exp< — - 



A (poi + P02Y 



pIi 



P02 



R' 



exp-^ 



which becomes 



TT 



exp 



(A+77)Poc + 



;pIs 



-^(A + i)(2po)^i?^ 



(38) 



(39) 



for the asymptotic singlet state. This is equivalent to the EFP of a spin-less 
Calogero-Sutherland system with coupling constant A' = X + 1/2, see (|43|) . 
This is consistent with the phase-space picture provided in in which it is 
explained that for a singlet state each particle occupies an area of 7r(A -I- 1/2) 
due to the exclusion statistics, while it would occupy an area 7r(A -I- 1) if it were 
alone. Therefore, in this context, the charge field can be thought of as describing 
a spin-less Calogero system with coupling constant A' = A -|- 1/2. 



5.1. Free fermions with 

Setting the coupHng parameter A = corresponds to non-interacting (free) 
fermions with spins and this reduces (|36|) to 



pfreefermions(^) = expj-l [tt (poc - Pe) i?]^ - J k (pos - P.) 



exp< [tt (poi - pi)Rf - ^ k {P02 -P2)Rf 



(40) 
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This result is the same as the one obtained in [17|. The EFP is then 

pfree fermions(^) = exp|-i {npo^Rf ^ (^Po2i?)^| , (41) 

which agrees with the results obtained in the context of Random Matrix Theory 
2^, 24 1, where the subleading corrections were also found. 

5.2. Spin-less Calogero- Sutherland model 

Of course, the prime check to our formula for the DFP/EFP of the sCM is 
to take its spin-less Hmit p2 = po2 ~ 0, which gives 

Pf;?-l^^^(i?) = eJ-^^n^po^-p^fRA, (42) 



(43) 



in perfect agreement with [17|, |2J| for a spin-less Calogero-Sutherland system 
with coupHng A' = A -I- 1. 

5.3. Probability of Formation of Ferromagnetic Strings 

If we require the minority spin particles to completely empty the region 
—R < X < R at T — 0, we are left only with the majority spin and we created 
a (partially) polarized state. We can refer to this case as the Probability of 
Formation of Partially Ferromagnetic Strings (PFPFS) [ii,E3|- Setting p2 = 
in ((36l) and leaving pi finite we have 

PpFPFs(-R) = exp|-y X{poi- pi+ po2f + ipoi- Plf + Pm -R^l • (44) 

The above is the probability of formation of ferromagnetic strings accompanied 
by a partial depletion of particles, since in the region of depletion we have pc = 
pi . We can impose that the average density of particles is constant everywhere 
by setting pi = poi + po2 = Poc while still requiring all particles in the region 
—R < X < R at T — to he completely polarized (maximal magnetization: 
PFFS) 

PpFFsiR) = expj- [7rpo2i?]'| • (45) 

Note that l(45l) is independent of A and exactly corresponds to the Emptiness 
Formation Probability ofaA' = A-|-l = 2 spin-less Calogero model with back- 
ground density given by po2 l[43| . Interestingly the same result l(52|) will be 
derived in the next sections as the EFP of minority spins, i.e. p2 = 0, in the 
Haldane-Shastry model (|46l) . This is just another aspect of the well-known re- 
lation between spin-Calogero, Haldane-Shastry and A' = 2 spin-less Calogero 
models [23, lisll as it will be shown in the next section. 
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5.4- The freezing limit 

If we take the A ^ oo limit in the spin-Calogero model the charge 
dynamics freezes (the particles become pinned to a lattice) and only the spin 
dynamics survives. This freezing limit was shown bv Polychronakos jH] to be 
equivalent to the Haldane-Shastry model (HSM) (Hlil: 



^^HSM — 2— -J ^ 



Sj ■ Si 



^ sin^ ^ U - I) 



(46) 



an integrable Heisenberg chain with long range interaction. In [3| the freezing 
limit was studied through a systematic expansion of the hydrodynamic fields in 
inverse powers of /i = A + 1/2: 



Pc Vc, Ps,Vs 



(47) 



With an additional rescaling of time t fj,t, the equations of motion were 
separated order by order in powers of p. 

It was shown that the charge sector is frozen, in that charge dynamics ap- 
pears only at orders 0(^~^) and higher, while the spin sector already has non- 
trivial dynamics at order 0(1). This dynamics is the same as the one derived 
independently for the HSM, over a background density of particles poc = N/L. 
As a consequence of charge freezing, this background density is kept fixed and 
constant up to order 0(1) and fluctuation are suppressed as Therefore, 
as fl — X + 1/2 oo, charge conservation is imposed dynamically everywhere, 
including in the region of depletion: 



-(0) ~(0) I ~(0) 

Pc ' ^ Pi +P2 = P01+ Pa2 = 



POc 



(48) 



The above l(48|l along with the usual depletion boundary conditions |(32 
reduces ll36l) to 



expi -y 



^(po.-pi"')%-p(^)+0(M-^) 



(49) 



exp 



'( 



Poi 



pf^]R 



= cxp 



TT (P02 - 4"') R 



where corrections for a flnite A are of the order 1 /p. Eq. (|49|) coincides with (|45|) , 
where condition (|48|) was imposed as a boundary condition, and not dynamically 
from the equations of motion. 



5.5. Haldane-Shastry model 

The hydrodynamic description of the HSM (|46| was constructed in (l9| . 
resulting in the following Hamiltonian for the minority spins (remember that 
the HSM is a lattice model and therefore there is no charge dynamics) 



HSM 



dx 



1 2 2 2 S 

2 P2 W2 + 3 TT P2 



(50) 
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which is the hydrodynamic Hamiltonian for a spin-less Calogero system with 
coupling constant A' = 2, see (fT3|) . It is in fact known that the spectrum of 
the HSM is equivalent to that of a spin-less Calogero-Sutherland model with 
exclusion parameter A' = 2, but with a high degeneracy due to the underlying 



Yangian symmetry [27 1 



The connection between the HSM l(50|) and the sCM in the freezing limit is 
to express the minority spin fields in (|50|) in terms of spin fields fl9| : 



POc - Ps „ ,_^s 

P2 = 2 ' = -2ws . (51) 

It is straightforward, using l(42|) with A' — 2, to see that the PFPFS for the 
HSM model is exactly l(49|) : 

P^I^FsiR) = expj- k (P02 - P2) . (52) 

The equivalence between (|45| , ijlQ]) and l(52|) is a strong check of the consistency 
of our methods and shows from a novel perspective the well-known relations 
between the sCM in the large A limit, the HSM and the spin-less Calogero- 
Sutherland model with A' = 2. 



6. Spin Depletion Probability 

So far, we considered a DFP problem specified by the boundary conditions 
l(32|) . i.e. by fixing the density of both species of particles on the segment 
r = 0, —R < X < R. However, our formalism allows for a more general and 
natural question. We can, for instance, demand a given magnetization (i.e. spin 
density) on the segment, without constraining the charge sector, i.e. imposing 
the boundary condition: 

PsiT = 0;-R<x<R)=ps , (53) 

instead of ((32l) . 

From ll34l) we have 



Ps ~Tr~~r:^'^ ~ ^^2] , (54) 

7r(A + i ) 

substituting in (|27|28|29p we find that ^ is satisfied if 

XaTja - XbTIb = XaPOa - hpob - (A + l)ps ■ (55) 

This equation leaves undetermined a complex constant ^ = + i^2 : for later 
convenience we parametrize the solution as 

XaVa = XaPOa - -j^Ps ^ ( + 2 ) ^ 
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^ + o (POc -0 + 77 (POs - Ps) , 



Ab?7b = Abpob + ( A + ^ ) (ps - 



(56) 



We have constructed the solution that realizes a constant spin density in the 
depletion region, while leaving the densities for the individual species free to 
vary. Please note that a finite imaginary part of ^ is necessary to have dxPiai'^ — 
0; -i? < X < i?) ^ 0. 

We can now substitute (l56l) in (ISTI) to find: 









{4 




-i) 







1 



(POs 



Ps) 



w 



(57) 

This probability depends on two, yet undetermined, parameters: ^1,2- If we 
choose ^1 = pc and ^2 = we recover exactly (|36l) . as we expected. This would 
correspond to forcing a given charge density at the depletion region, together 
with (pl) . 

We also note that the configuration that maximizes the probability l(57|) is 
given by ^ = poc: 



PfSUR) = PsDp{R; e = Poc) = exp^ -— {pos - PsY 



(58) 



One cannot help but noticing the similarity between l(58|) and (|49| or l(52|) . 

Since the parametrization we choose in ((56l) allowed us to express the prob- 
ability l(57|) as a Gaussian for ^1,2, we would get the same result by performing 
an integral over the free parameters: 



00 J —00 



(59) 



where the prefactor coming from the Gaussian integration is beyond our accu- 
racy anyway (note, however, that it does not depend on the charge density). 
The integration over ^ corresponds to summing over all configurations of the 
form l(27l [28112911. 

The probability of realizing the magnetization set by ^53^ is given by a sum 
over all configurations that satisfy the given boundary conditions. To perform 
this sum correctly, we would need to consider all possible charge density profiles 
Pcix) at the depletion region and therefore consider more general solutions than 
I . These general solutions are of the form 



AaPOa — -^Ps 



7TPs+^ A+- C(^), 



(60) 



AbPofc + ( A + ^ I p. 



i?2 



+ 7r A + 



-Ps+Ciz)] 
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where ^(z) is an analytic function such that Re^(x) — Pc{x) and ^(z oo) 0. 
The sum over all configurations satisfying l(53|) can be formulated as a functional 
integral over all functions £,{z). However, it is easy to convince oneself that the 
configurations that minimize the action are of the form l(28l[29l) . with ^ — poc- 



7. Charge Depletion Probability 



The last problem we will address is conjugated to the one considered in the 
previous section, i.e. the probability of reaHzing a given depletion of the charge 



Pc{t = 0;-R< X < R) ^ Pc , 
without constraining the spin density. Using (|34l) we have 

Pc = Re 



which means that (|27l28l29p fulfill ^ if 

Va+Vb = POa + POb - Pc 



(61) 



(62) 



(63) 



Once again, we are left with the freedom of introducing a complex number 
^ = ^1 + i^2 to parametrize the solution: 



2A+ 1 



1 



Vb 



POa 



2{X + lY'^ 2(A + 1) 



2A+ 1 
2(A+1) 



(POc - Pc) + 
1 



1 



2(A + 1) 



2(A+1) 

Inserting this into ((3T1) we obtain: 



2(A + 1)^^- 2(A + 1) 

(pOc - Pc- POs + C) 



(64) 



Pci3p(-R;0 =cxp - 



1 



1 



(65) 

Setting ^1 = ps and ^2=0 correctly reproduces l(36|) . while the maximal prob- 
ability is achieved for ^ — pos'. 



PCDP{R) - PCDP{R: e - PO.) = CXp 



+ ^ ) (POc - Pcf 



R" 



(66) 



As before, since (|65l) is Gaussian in ^1^2, we would obtain the same result by 
integrating over these variables 



/oo ncx: 
00 J —00 



(67) 
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where we neglected the coefficient coming from the Gaussian integration because 
its beyond the accuracy of our methodology. This integration corresponds to 
summing over all configurations given by (|28l [29l) . Again, we notice a striking 
similarity between l(66l [67|) and ((36l) . We will comment in the next section on 
how to interpret these results. 



8. Discussion of the results 



Eq. l(36|) looks like the product of the two independent depletion probabilities 
for the spin and charge sector. This interpretation is supported by the results 
of the two previous sections, see l|58p and l|66p . but it is quite surprising in a 
sense. In fact, it would indicate a sort of an effective spin-charge separation, as 
if the spin and charge degrees of freedom could be depleted independently. This 
is contrary to intuition, since spin-charge separation is realized only for low- 
energy excitations close to the Fermi points, while the EFP involves degrees of 
freedom deep within the Fermi sea (and requires a full non-linear hydrodynamic 
description beyond the usual bosonization approach). However, it seems that 
from a EFP perspective spin-charge separation survives beyond the linearization 
of the spectrum, at least at leading order for the sCM. 

In fact, quite surprisingly, for the Calogero-type interaction (as well as for 
free fermions), the DFP result obtained for small depletions using a linearized 
hydrodynamics (conventional bosonization) can be extended up to a complete 
emptiness and remain quantitatively correct This is due to the fact that 
the gradientless hydrodynamic for this interaction is purely cubic, see (fTTj) . 
This fact has two important consequences: the first one is that the equations 
of motion can be written as Riemann-Hopf equations (|25l) , which are trivially 
integrable with the implicit solution given by l(26|) . This is important to connect 
the boundary conditions at infinity with those due to the DFP. The second fact 
is connected to the form of the parameters u and k of the linearized theory, 
which in Hamiltonian formalism can be written in general as 



H = dx 



U „n UK,, .21 . , 

+ — (V0)' , 68 

where n(a::) and 0(x) are conjugated fields. In terms of hydrodynamic variables 
l(20|) they are 

v{x) = n{x) , p{x) = po + V0(a;) , (69) 

where po is the background value over which we are linearizing the theory. Note 
that K = where K is the conventional Luttinger parameter [sot . 

For Calogero-type models, the sound velocity u depends Hnearly on the 
density, while the interaction parameter n does not depend on the point around 
which we are Hnearizing. In appendix[B]we show that the sound velocity can be 
rescaled out of the DFP calculation (at zero temperature) and k is the relevant 
factor encoding the interaction, which determines the coefficient of the Gaussian 
behavior of the DFP. All these peculiarities of the Calogero interaction conspire 
in a way that extending the small depletion result to higher depletion is "trivial" 
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and, in fact, gives the correct result. Let us remark in this respect, that any non- 
Hnear theory can be seen as the integration of successive Hnear approximation, 
where the coefficients are adjusted at each point. In this Hght and from what 
we pointed out above, it is clear that the simpHcity of the Calogero interaction 
allows a simple integration of successive linear theories for the DFP calculation 
and this is the reason for which the linearized result can be trivially extended 
from a small DFP to a complete EFF. 

In appendix[B]we calculate the DFP in the Hnearized approximation (guided 
by 0). The result is lfT05l) 



S 



linear 
DFP 



(70) 



where we used rj = fj = pQ — p. 

If we substitute lfT2l) in (fTTI) we can write the hydrodynamic Hamiltonian in 
terms of spin and charge fields as 



H = 



dx ■ 



1 



-,PcVc 



TT 

~6 



A + 



A + 2 ) - Ap, 



PsVcVs 



TT- / , 1 \ ^2 ^ X „3 



(71) 



- A + -lp.p^--Ap^ 



By linearizing this Hamiltonian, i.e. by expanding the fields as Pc,s — Po;c,s + 
Spc,s and looking at the coefficients in front of the quadratic part, we find the 
following parameters: 



TT ( A + i ) poc , 



UciPOc) 



Usipos) = TT A + - poc - Ti'A/Co 



KcipOc) = TT A + 



1 



KsiPOs 



TT 

2 



(72) 
(73) 



We see that substituting these values in ((TOl) correctly reproduce i(58t [66|) and 
therefore l|36p as well. This means that not only the linearized theory is sufficient 
to calculate the correct coefficients of the EFP, but also that the spin-charge 
separation survives as if the linear theory was valid for high depletions as well. 

To conclude, we can suggest a simple physical interpretation of l(70|) . In 
a Calogero-Sutherland system, the interaction parameter k — ttX' has a simple 
semiclassical interpretation in terms of the phase-space area occupied by a single 
particle, see 31 1 and We can then see that (|70| represent a volume in the 
x — T — k space: the phase-space area at a given r is of the order of k{po — p)R, 
see, for instance, (|24|) . This has to be multipHed by the number of particles 
involved in the depletion over time, which is of the order (po — p)R- 



9. Conclusions 

We calculated the Emptiness and Depletion Formation Probability for the 
spin Calogero-Model ^ and for the Haldane-Shastry Model {Ml). The EFP 
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is one of the fundamental correlators in the theory of integrable models and, 
despite its being non-local, is considered to be one of the simplest. Nonetheless 
its asymptotic behavior is known only for a few systems and in this paper we 
calculated it for the sCM and the HSM for the first time, at the leading order. 

The DFP is a natural generalization of the EFP in the hydrodynamics for- 
malism we employ. By calculating the DFP in its most generality (|36|) . we 
can achieve the different EFPs by taking its appropriate limit. The calculation 
is done in an instanton picture, where the DFP is viewed as the probability 
of formation of a rare fluctuation in imaginary time that reaHzes the required 
depletion at a given moment. 

The long distance asymptotics in 1-D models are normally calculated in a 
field theory approach using bosonization. However, as this approach is valid 
only for low-energy excitations close to the Fermi points where the linearization 
of the spectrum is a reasonable approximation, it is not sufficient for the EFP, 
which involves degrees of freedom deep in the Fermi sea. For this reason we 
used a non-linear version of bosonization, i.e. the hydrodynamic description 
developed in (l9| . 

All our formulae show a characteristic Gaussian behavior as a function of 
the depletion radius R. This is to be expected for a gapless one-dimensional 
system, as it was first argued in This is because in the asymptotic limit we 
consider, R is the biggest length scale in the system and therefore the instanton 
configuration will have a characteristic area of i?^, where the second power comes 
from the dimensionality of the space-time. In [3;[l3| it was also shown that for 
small depletion, the linearized bosonization approach is sufficient to calculate 
the DFP, while in general it deviates from the correct results for progressively 
bigger depletion and, eventually, emptiness. 

However, the Calogero-Sutherland kind of models (as well as non-interacting 
fermions) are special and the linearized result happen to coincide with the cor- 
rect, non-linear one. We argued on the origin of this observation in the previous 
section. Moreover, we noticed that l(36|) and the analysis of section [6] and [7] 
indicates that, from a EFP perspective, spin-charge separation seems to sur- 
vive beyond the linear approximation, in disagreement with what one naively 
would expect. This resurgence of linear results in a non-linear problem is a very 
surprising result, peculiar of the sCM. 

The coefficients in front of R^ are novel of this work. In section [8] we inter- 
preted them from a bosonization point of view and via a simple semiclassical 
argument and throughout the paper we have checked them against known results 
in certain limits where possible. In particular, we showed agreement with the 
free fermionic limit l|40p and the spin-less Calogero-Sutherland model l(42|) . For 
both of these models, the EFP has a particular interest coming from Random 
Matrix Theory, as it is known that for certain rational values of the coupling 
parameter A the CSM describes the RMT ensembles. It would be interesting if 
the sCM would also have an interpretation in terms of some generalized random 
matrix model, but we are not aware of such connection yet. 

In section 15.41 we used the fact that the Haldane-Shastry model can be 
achieved as the freezing limit (A oo) of the sCM to calculate the Probability of 
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Formation of (Partially) Ferromagnetic Strings in the HSM. In section ISTSl the 
same quantity was derived independently from the hydrodynamic description 
of the HSM. Section [5^ and [531 highlight the correspondence between large- A 
sCM, HSM and A' = A + 1 = 2 spin-less Calogero model from a EFP/DFP 
perspective. 
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A. Action for the DFP solution 

In this appendix, we will revise the calculation presented in and adapt 
it to our case. We want to find the value of the hydrodynamic action calculated 
on a given solution satisfying the DFP boundary conditions. 

The gradientless hydrodynamic action in imaginary time t = it can, in 
general, be written as 

S = ld\C[v,p]^ |dxdTp|y + , (74) 

where 

<P) - (75) 
is the internal energy per particle of a Calogero system and 

P^d,[pe{p)L^^^^^pl (76) 



is the chemical potential. The action l(74|) has to be supplemented with the 
continuity equation 

drp + (pv) = , (77) 

which can be considered as a constraint relating the two conjugated fields p and 
V. This constraint can be resolved by introducing the displacement field (j){x, r): 

p = Po + dx<p , j = pv = -drcj) ■ (78) 
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Physically, the displacement field counts the number of particles to the left 
of a point. We can use l(78|) to write the Lagrangian as a functional of 4>: 
C[p,v] = Its variation then gives the Euler equation for the fluid, which 

can be written more simply as 



drV + vdxV = dxdp [pe{p)\ ^ \ n pdxP ■ 



(79) 



For the particular choice of internal energy l(75l) . corresponding to the Calogero- 
Sutherland interaction or exclusion statistics, the Euler equation and the conti- 
nuity equation can be combined into a single complex Riemann-Hopf equation: 

drk — ikdxk = , fc(T, x) = A7rp(T, x) + w{t, x) . (80) 

This equation has a simple, implicit, solution of the form 

k = F{x + ikT) . (81) 

In the body of the paper, we argued that a solution satisfying the DFP 
boundary conditions is of the form 



F{z) = F{z; po,r]) = Xirpo + Xttt] 



1 



(82) 



Here, po is the background (equilibrium) density at infinity (where moreover 
V = 0), and rj is a, possibly complex, constant specifying the DFP. 

To calculate the Depletion Formation Probability, we need to compare the 
action ([74|) calculated on the solution (|81l82p to the action of an equilibrium 
configuration: 



dx dr <j p— + pe{p) - poe(po) - p{p - Po) 



(83) 



To take advantage of the fact that ([81]) is a solution of the equations of motion, 
we first take the variation of (|83|) with respect to the parameters of the solution. 
In this way, we will reduce a two-dimensional integration to a contour integral 
over the boundaries, since the bulk terms are proportional to the Euler-Lagrange 
equations and vanish: 



d (5 - Sq) = dp„ {S - 5o) dpo + 9^ (5 - Sq) d-q + df^{S - So) dr] 
We have: 



(84) 



dr, (S - So) 



d X 



d X 



-vdr4>n 



dp (pe) + p 



dx(l>v 



dr [V C/Jri] - dx 

drV + vdxV — dxdp (pe) 



Y - 5p (pe) 



[v dx 



dp (pe) 



dt 



(85) 
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where 4>ri = dr^4>- 

The boundaries over which the contour integral is taken are, by Stoke's 
theorem, the points where the integrand has a discontinuity. It is easy to check 
that this contour comprises only two paths: one at infinity (Co = + i/crj = 
cx)}) and one around the branch cut of (|82|) . which we take along the real axis 
(Ci = {r = 0±,-i? < X < R}). 

From l|81|82p . at infinity we have 



p(ko| - 


-f oo) 




v{\zq\ - 


-> oo) 




0(ko| - 


-> oo) 


4 Uo 



(86) 
(87) 
(88) 



and we see that the integrand in l(85|) along Co vanishes too fast and the contour 
integral gives no contribution. 

Close to the cut on the real axis we have 



^0 ^0 

V__ J_ 



p{r = 0^ 



-R < X < R) = Pa 



=F 1 



2 y/R^ - 



v{t = 0* 
<I>{t = 0± 
Therefore 

dr, {S - So) = 

Similarly, we have 

df^ {S - So) = 



-R<x < R) 
-R<x <R) 



lAvr ^ ATT 



2 



— ^ — X ± 1 — - — 



(89) 
(90) 
(91) 



R 

[v{x,Q+)(f>^{x,Q+) - v{x,Q-)cl),^{x,Q-)] dx 

-R 

Xn^R^ _ 



-V 



[v (f>fj] dx 



(92) 



dp (pe) + U-r; 



dt 



Xn^R^ 



-V 



(93) 



The derivative with respect to po is a bit more complicated as it involves 
more terms. After a bit of algebra and an additional integration by parts we 
obtain: 



dpo {S -So) 



dx 



dp (pe) + p 



+ x) + (5po^)0 



dtj (94) 
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where 4>p^^ = dp^ 4>. Substituting the behaviors (|86ll88p and (|89ll9ip , the integrals 
around the two contours gives equal but opposite results (±^A7r^(?7 + ffjE?), 
which means 

dp, (5 - 5o) = , (95) 

as one could have expected. 

We can now integrate ^ using (|92l93|95p to find 

5dfp -5-5o - iATr^TyTyi?^ . (96) 



B. Linearized Hydrodynamics and DFP 

Let us linearize the hydrodynamic equations, by expanding the theory and 
retaining only the quadratic part of the Lagrangian. This linearized hydro- 
dynamics is usually referred in the literature on one-dimensional models as 
hosonization. 

From the previous section, we have that the gradientless hydrodynamic La- 
grangian for a one-component system is 



2p 



(97) 



where j = pv. We expand the fields around a background value and we 
parametrize the fluctuations around this background through the displacement 
field (j) as in l(78|) . so that the constraint l(77|) is automatically satisfled. Keeping 
terms up to the quadratic order we have 



1 

2^ 



1 



(dr^r + T^d'p {P^)p=p, (d.^by + po e{po) - p 



-{OA) + — {d., 



(98) 



where in the last line we introduce the standard parameters of bosonization: 
the interaction parameter k — k{po) and the sound velocity u — u{po). The 
displacement field evolves according to a linear wave equation: 



9^0 + u^dl(l> = 



(99) 



The linearized treatment is valid for small fiuctuations around the back- 
ground po, i.e. as long as the gradients of (f) are small and only low energy 
excitations are involved. For this reason, it is not possible to calculate the EFF 
through standard bosonization, but we can consider a DFP with very small 
depletion. 

It is simple to see that the solution that satisfies the DFP boundary 
conditions is of the form: 



Re 



(100) 
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where zq = x + m{pQ)T. For this solution to be compatible with the linearized 
approximation we need \r]\/po ^ 1. 

It is easy to calculate the DFP by evaluating the linearized action 



S-Sn 



I dxdr {±{d,cj,f + ^ } (101) 



on the solution (|100p . One immediately observes that, at zero temperature, this 
action does not depend on the sound velocity, as we can rescale the time as 
y = ut: 

S-So^'^j dxdy [{dy<j,f + (9,0)'} . (102) 
We can further rescale the lengths by R and substituting (|100p in (|102p we get 

5 - 5o ~ I Wi?' j didy \^\x + ij/)|' , (103) 

where 

0(z) = \/z2 -I- z . (104) 

Now all the physical parameters have been explicitly extracted and one has just 
to perform an integral that contributes only with a numerical factor. The result 
is 

S -So-^nriT^R^ . (105) 

Since for a Calogero-Sutherland system k — Att, we notice that (|105p exactly 
coincide with (|96l) . This is quite surprising, since, as we argued above, the 
Hnearized result should be trusted to be approximately correct only for small 
depletions. However, the Calogero kind of interaction is very special and we 
can extend l|105p to higher depletion, without loosing accuracy. In section [8] 
we discuss the meaning of this observation. Let us remark that this result is 
specific for the EFP and it is not to say that for Calogero-Sutherland systems 
the effects of non-linearity are in general not important. For instance, effects of 
non-Hnear spin-charge interactions were observed and discussed in [19l |. 

This DFP calculation can also be performed using the line integral technique 
explained in the previous section. In this case, the variation of the action (|10ip 
gives simply: 

dri{S - So) ~ ^ { M ^^''^^ 4>rjdx + Ku {dx4>) (j)r, dr 

= K(j)l^{dy<j))(l),^dx + {d^<j>)(f>,jdyY (106) 
One can then proceed as we showed in the previous section to easily recover 
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